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ABSTRACT 


INTRODUCTION 


A probability-based analysis method for 
predicting buckling loads of compression-loaded 
laminated-composite shells is presented, and its 
potential as a basis for a new shell-stability 
design criterion is demonstrated and discussed. 
In particular, a database containing information 
about specimen geometry, material properties, 
and measured initial geometric imperfections for 
a selected group of laminated-composite 
cylindrical shells is used to calculate new 
buckling-load "knockdown factors.” These 
knockdown factors are shown to be substantially 
improved, and hence much less conservative 
than the corresponding deterministic knockdown 
factors that are presently used by industry. The 
probability integral associated with the analysis 
is evaluated by using two methods; that is, by 
using the exact Monte Carlo method and by 
using an approximate First-Order Second- 
Moment method. A comparison of the results 
from these two methods indicates that the First- 
Order Second-Moment method yields results 
that are conservative for the shells considered. 
Furthermore, the results show that the improved, 
reliability-based “knockdown” factor presented 
always yields a safe estimate of the buckling 
load for the shells examined. 
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The difficulty in predicting the behavior of 
compression-loaded thin-walled cylindrical shells 
continues to concern design engineers. Thanks 
to extensive research programs conducted in the 
late sixties and early seventies of the past 
century and the contributions of many eminent 
scientists, it is now known that buckling strength 
can be affected significantly by the uncertainties 
in the definition of loads, material properties, 
geometric variables, edge support conditions, 
engineering models, and the accuracy of the 
analysis tools used in the design phase. 

Presently, the NASA design criteria monographs, 
e.g. NASA SP-8007 1 , from the late 1960’s are 
used, to a wide extent, in industry to design 
buckling-critical shells. The data in these 
monographs account for design uncertainties by 
using a “lump-sum” safety factor that is 
determined from experimental data. This safety 
factor is often referred to by designers as the 
“empirical knockdown factor y Structural 

design verification is achieved by applying the 
worst-case loading to the structure and testing to 
failure. However, it has been found that using 
this approach to design buckling-critical shells 
often leads to overly conservative designs. The 
reason for this deficiency is primarily due to the 
fact that the empirical knockdown factor is based 
on a lower-bound design curve derived from a 
large amount of experimental data for which the 
test pedigree and other things like initial 
geometric imperfections and boundary 
conditions are not well defined. In addition, 
these knockdown factors do not include 
information for shell structures made from 
advanced composite materials. Thus, new, 
improved knockdown factor are expected to yield 
significant weight savings and improved 
performance in the design of buckling-critical 
shells. Further, the new knockdown factors 
presented herein will, for the first time, provide 
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design information for compression-loaded 
composite laminated shells. One approach that 
appears to have potential for realizing these 
improvements is a reliability-based approach. 

Following the pioneering work of Freudenthal 2 on 
structural safety and reliability, around 1970, the 
use of probabilistic methods to assess the safety 
of structures was proposed. 34 Recent papers, 
e.g. Refs. 5-7, have proposed a new reliability- 
based design procedure for buckling-critical 
isotropic shells. However, similar results for 
compression-loaded laminated composite shells 
do not exist in the literature. 

The focus of the present paper is on the 
extension of the probabilistic approach to 
buckling of compression-loaded laminated 
composite shells. In particular, the goal is to 
formulate a new and improved scientific 
knockdown factor that is based on statistical 
information obtained from measurements made 
on families of composite shells. As a minimum, 
the measurements include initial geometric 
imperfections, but can also include other critical 
details such as shell-wall thickness variations, 
material property variations, non-uniform load 
introduction, and boundary condition effects. 
With this approach, it is believed that a 
scientific, reliability-based knockdown factor will 
result in a better engineered, better designed, 
and safer structure by quantifying, and supplying 
understanding of, problem uncertainties such as 
initial imperfections and their interaction with the 
design variables. 

Thus, the objective of the present paper is to 
demonstrate a reliability-based buckling-analysis 
approach that can be used as the basis of a new 
preliminary design criterion for buckling-critical 
shells, particularly, for shells made of advanced 
laminated-composite materials. To accomplish 
this objective, a probability-based analysis 
method for buckling of compression-loaded 
cylindrical shells is presented first. In this 
method, a database that contains information 
about test specimen geometry, material 
properties, and measured initial imperfections for 
a selected group of laminated-composite 
cylindrical shells is used to calculate reliability- 
based buckling-load knockdown factors. In 
addition, the probability integral associated with 
the analysis is evaluated by using an 
approximate, less computationally intensive, 


First-Order Second-Moment (FOSM) method, 
and these results are compared to 
corresponding results obtained by using the 
exact, computationally intensive Monte Carlo 
method. Then, the results of the reliability-based 
analysis are compared with results based on 
NASA design guidelines for buckling of 
compression-loaded isotropic cylindrical shells. 


RELIABILITY-BASED BUCKLING ANALYSIS 

The buckling problem for axially compressed 
cylinders can best be formulated in terms of a 
response (or limit state) function such as 

g(X) = \ s (X)-X (1) 

where A is a suitable normalized loading 
parameter (e.g., A = P/P C ), A s is the random 

collapse load of the shell, and the vector X 
represents the random variables of the problem. 
The components of the random vector, X h may 
be Fourier coefficients of the initial geometric 
imperfections and other parameters that quantify 
the uncertainties in the specified boundary 
conditions, the constitutive equation used to 
describe the nonlinear material behavior, the 
shell-wall thickness distribution, residual 
stresses, etc. Notice that the evaluation of the 
response function generally involves the solution 
of a complicated nonlinear structural analysis 
problem that typically requires a detailed, and 
possibly large, finite-element model. However, 
with today’s computational resources, a complex 
nonlinear analysis in itself does not pose any 
Insurmountable difficulties. The response 
function g(X) = 0 separates the variable space 
into a “safe region,” where g(X)> 0, and a 
“failure region,” where g(X)< 0. The reliability 
R( A) is calculated from the probability of failure 

PA A) by 

AW=1-^(A) (2) 

where 

PfW = Prob {g(X) < 0} = L f f*(x)dx (3) 
dm o 

and fy(x) is the joint probability density function 
of the random variables involved. The credibility 
of this approach depends on two factors; the 
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accuracy of the mechanical model used to 
calculate the limit state function and the 
accuracy of the probabilistic techniques used to 
evaluate the multi-dimensional probability 
integral. 

The First-Order Second-Moment (FOSM) 
Method is an approximate method for evaluating 
multi-dimensional probability integrals that is not 
as computationally demanding as exact methods 
such as the Monte Carlo method. The FOSM 
method involves linearizing the response 
function such that g(X) ~ Z(X) at the mean point 
and requires knowledge of the distribution of the 
random vector X . By assuming that both Z[X) 
and X are normally distributed, then 

/^A) = 1-Prob(Z< 0) 

f (4) 

= 1- J f z (f)dt K ’ 


where the probability density function f z {t) is 
given by 


= 


a z ^ 2 k 


exp 


1 fW 


(5) 


and a=E(Z), <r z = A /var(Z). It has been shown 
in earlier works (e.g. Refs. 5 and 6) that with this 
approach, the solution of the reliability-based 
buckling problem defined by Eq. (3) can be 
reduced to a series of /7+1 deterministic 
buckling-load calculations, where n is the 
number of random variables used. 


amplitudes, respectively; and the axial and 
circumferential wave numbers /;A,and^ are 
integers. 

From the work of Koiter 8 and Budiansky and 
Hutchinson 9 , it is known that the worst 
imperfection shapes (the ones that give the 
biggest decrease in the buckling load for a given 
amplitude) are those which are similar to the 
critical bifurcation buckling mode-shapes. Thus, 
predicted bifurcation buckling mode shapes are 
used in the present analysis to form a worst-case 
imperfection shape. In order to apply the First- 
Order Second-Moment method, the amplitudes 
of the axisymmetric and asymmetric components 
in the imperfection model given in Eq. (6) are set 
equal to the root-mean-square values of the 
corresponding measured initial imperfections. 
The initial imperfection is modeled using the 
following alternate double Fourier representation 


/»i 

W(x,y) = t2 j s4 0/ cosm - 

/= i L 


(7) 


+'X X sin Ax - C k , cos i - + D k , sin l ^ 

k = 1 (=2 ^ ' 


Then, by definition, the root-mean-square value 
of the measured initial imperfection is 


2izR L 


0 0 

/?-) n ^ /?2 


(8) 


-T 2S4+XK4+4,) 

/=1 k=\ (=3 


NUMERICAL RESULTS 


To include the effects of initial geometric 
imperfections in the analytical model in a simple 
manner, for the purpose of demonstrating the 
reliability-based approach, the following two- 
mode initial imperfection model that consists of 
an axisymmetric and an asymmetric component 
was employed: 

W — x — x v 

— = ^cos/7r j + ^sinAw— cos£-^ (6) 

where W is the out-plane imperfection; R, L, 
and t denote the shell radius, length, and 
thickness respectively; f, and are the 
axisymmetric and asymmetric imperfection 


Thus, 


f 


h , n 2 


2 ,=i 
= ^axi + ^asy^- 


k = 1 f-3 


(9) 


The measured initial shape of one of the 
composite shells used (AWCYL-11-1) is shown 
in Fig. 1 . As can be seen, the initial imperfections 
are dominated by the £= 2 mode, which 
represents an ovalization of the cylinder. 
However, the eigenvalue from a linear bifurcation 
analysis corresponding to the asymmetric mode 
(k= \ t=2) is equal to 4.106, which is about 10 
times bigger than the lowest eigenvalue 
X c =0.398. Based on this fact, it was decided to 
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neglect the contribution of the ovalization (the 
( = 2) mode when calculating the root-mean- 
square values of the measured initial 
imperfections. Otherwise, the resulting numerical 
model used for the buckling load calculations will 
predict unrealistically low buckling loads because 
the amplitude of the imperfection would be 
artificially inflated due to the large magnitude of 
the ovalization (the t = 2) mode. 

In this paper, the values |, = -A axi and f 2 = & aS y 

denote the axisymmetric and asymmetric 
imperfection amplitudes, respectively. A listing of 
the imperfection amplitudes obtained for nine 
composite laminated shells fabricated and tested 
at NASA Langley Research Center are given in 
Table 1. The buckling load calculations were 
carried out with the code TWOMOD 10 . This code 
is based on a Galerkin-type solution of the 
nonlinear Donnell-type equations for imperfect 
anisotropic cylinders. The theoretical buckling 
load is defined as the value of the 
nondimensional loading parameter X = PI P d at 
the limit point of the prebuckling state. For the 
statistical calculations, the data associated with 
laminated composite shells are used. A listing of 
the geometric and material properties of the nine 
NASA shells tested is given in Table 2. Typical 
results of the data reduction for the measured 
initial imperfections of theses shells are 
tabulated in Refs. 11 and 12. The mean values 
and the variance-covariance matrices of the 
random variables Xj were evaluated by using 
the following ensemble averages for a sample of 
the experimentally measured initial imperfections 


For the buckling load calculations, the shell 
AWCYL-11-1 was used. The imperfection model 
of Eq. (6) was specialized to the most critical 
long-wave imperfection combination with the 
axisymmetric mode /= 2 and the asymmetric 
mode A'=1and£=6. In particular, the strong 
coupling condition i=2k is enforced for this set 
of imperfection mode-shapes. This condition for 
strong coupling between an asymmetric and an 
axisymmetric imperfection shape was first 
identified by Koiter and has been shown to 
cause a significant reduction in the buckling load 
of a compression-loaded cylinder. The 
corresponding eigenvalues are X ci = 5.150 and 
X cM =0.401. In the probabilistic analysis, it is 
convenient to re-normalize these eigenvalues so 
that the critical (the lowest) re-normalized 
eigenvalue is equal to 1.0. Thus, 


Pci = 


0.398 


= 12.932 


and 


Pc, At - Q 


cM 


398 


= 1.007 


In order to apply the First-Order, Second- 
Moment (FOSM) Method, the mean buckling 
load has to be calculated first. Using the 
imperfection model of Eq. (6), with the mean 
values given in Eq. (12), yields E(A S ) = 0.896. 


Thus, 


E(Z) = E(A s )-p = 0.896 -p (14) 




777=1 


( 10 ) 


where p is the re-normalized load parameter 
defined as 


cov(X J ,X k ) = 


1 

M - 1 


M 


/ 77 = 1 


( 11 ) 


yielding 


E(X^) =-0.00334 
E(X 2 ) = -0.0498 


co\z(X v X 2 ) = 

’ 0.0000108 -0.000101 
-0.000101 0.00121 


(13) 


p= — (15) 

K 

and X™ = 0.398 is the critical (lowest) normalized 
buckling load of shell AWCYL-11-1 used for the 
numerical calculations. The variance Z is given 
by 


var(Z) = var(AJ 
2 2 

-II 

j-\ k = 1 


'3A. V 




'i ) 


3A< 

Wa 


cov(X v X 2 ) 


(16) 


The calculation of the derivatives 8A S I8%, is 
performed numerically by using the following 
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numerical differential formulas evaluated, at 
values of = E{Xj ) : 

3A, A s & + A§ t ,S 2 )-A s fo,&) (1?) 

3§i A £i 

8 &- ^ <8) 

For the increment of the random variable, 1% of 
the original mean value of the corresponding 
equivalent Fourier coefficient is used, so that 
A = 0.01 E{Xj) . The calculated derivatives are 

^- = 0.748 and ^- = -1.362 (19) 

"Si dg 2 

Next, using the sample variance-covariance 
matrix of Eq. (13), the variance of Z is calculated 
as 


knockdown factor y for isotropic shells, and the 
corresponding experimental buckling loads /| xp 
is presented in Table 3. As can be seen from 
the results, the reliability-based improved 
“knockdown” factor p a always yields a safe 
estimate of the load-carrying capacity of the 
shells tested if failure is caused by elastic 
instability. In addition, the results indicate that 
the new knockdown factors can be much less 
conservative than the corresponding traditional 
empirical knockdown factors for these shells 
based on the NASA design monograph for 
isotropic shells. For example, The traditional 
design approach for shell AWCYL-1-1 predicts a 
buckling load 36% lower than the experimentally 
measured buckling load. In contrast, the new 
probability-based knockdown factor predicts a 
buckling load that is only 7% less than the 
experimental buckling load. 

CONCLUDING REMARKS 


var(Z) = 0.00246 (20) 

Hence, the reliability index /? is 

F(Z) £-(A,)-p _ 0,896 -p 

-y/var (Z) °z 0.0496 

Finally, the reliability is calculated directly from 
Eq. (4); that is, 

Atp) = 1-P A (p) = 

1 ( 22 ) 

1-^z(0) = 2 + eff(P) = m 

where </>(/?) is the standard normal probability 
distribution function. This result is plotted in Fig. 
2. Notice that for a reliability of 0.99999, one 
obtains an improved, re-normalized “knockdown” 

factor of p a =X a l X m c = X a \ 0.398 = 0.68. 
Results from the FOSM are compared to similar 
results obtained from a Monte Carlo method in 
Fig. 2. The Monte Carlo method is based on 
5000 realizations or sample points as compared 
to the three analyses required to perform the 
FOSM method. These results indicate that the 
results obtained from the simplified FOSM 
method are conservative in the high reliability 
region the reliability curve, e.g., R(p) greater 
than 0.4. 

A comparison of the improved, scientific 
knockdown factor p a , the traditional empirical 


A probability-based analysis method for 
predicting buckling loads of compression-loaded 
composite cylinders has been presented. It has 
been demonstrated that a database containing 
information about the specimen geometry and 
material properties and measured initial 
imperfections for a group of specimens can be 
used successfully to calculate improved, less 
conservative, buckling load “knockdown factors,” 
as compared to the lower-bound knockdown 
factors currently used in industry design practice. 
A First-Order Second-Moment (FOSM) method 
and a Monte Carlo method were both used to 
evaluate the probability integral and a 
comparison of the results from both methods 
indicated that the FOSM method can be 
conservative for high values of reliability. 
Furthermore, the results show that the reliability- 
based “knockdown” factor p a always yields a 
safe estimate of the buckling load of the cylinder 
specimens if failure is caused by elastic 
instability. 

The probability-based analysis procedure 
presented in this study can be used to form the 
basis for a shell analysis and design approach 
that includes the effects of initial geometric 
imperfections on the buckling load of the shell. 
In particular, a probability-based analysis 
procedure can address some of the critical shell- 
buckling design criteria and design 
considerations for stability-critical shell structures 
without resorting to the traditional empirical shell 
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design approach that generally leads to overly 
conservative designs. 

It is anticipated that, for applications where the 
total weight of the structure is one of the critical 
parameters (e.g., aerospace structures), there 
will be a chance for definite improvement in the 
design process with the use of a new 
probabilistic design procedure. It is felt that the 
small added effort involved in systematically 
carrying out the required initial imperfection 
surveys will be fully justified by the overall cost- 
savings and by producing improved, less 
conservative and more reliable shell structures. 
However, in order to develop a validated 
probabilistic design procedure for buckling- 
sensitive structures, there is a need for additional 
systematic combined experimental and analytical 
or numerical results. Only then can the benefits 
be reaped from the many years of concentrated 
shell buckling research of the late 1960’s and 
early 1970’s. It is the authors’ opinion that the 
technology now exists for such an undertaking. 
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Table 1. Values of root-mean-square imperfection amplitudes. 



<?1 ~ ~&axi 

£2 “ A asy 

AWCYL-1-1 

-0.27460 -lO -3 

0.15395-10" 1 

AWCYL-2-1 

-0.33547 -lO" 2 

0.6891 5 -10 _1 

AWCYL-3-1 

-0.23708-1 O' 2 

0.30062-1 0 _1 

AWCYL-4-1 

-0.40 151-1 O’ 2 

0.37061-1 0' 1 

AWCYL-5-1 

-0.28099 -1CT 2 

0.321 79 1 0 _1 

AWCYL-1 1-1 

-0.22241 10 -2 

0.34885 -10 -1 

AWCYL-92-01 

-0.25832-1 O' 2 

0.381 90-1 0' 1 

AWCYL-92-02 

-0. 1 1 554 • 1 0 _1 

0.13259 

AWCYL-92-03 

-0.89050 -lO" 3 

0.58660 -10 -1 


Table 2. Geometric and material properties of the NASA composite shells. 6 


‘ 


R[in] 

t[in] 

tply pH] 

E u [/6/tf] 

E 22 [/t >l/rf] 

G^lblirl 1 ] 

Vi 2 

AWCYL-1-1 

[ ±45 / 0 / 90] 5 

7.999 

0.0399 

0.00499 

18.51 1-10 6 

1.64 • 10 6 

0.871 -10 6 

0.300 

AWCYL-2-1 

[±45/+45] 2s 

8.009 

0.0797 

0.00497 

18.578 10 6 

1.64- 10 6 

0.874 10 6 

0.300 

AWCYL-3-1 

[±45/0/90] 2i 

8.004 

0.0793 

0.00495 

18.671 10 6 

1.64 10 6 

0.878 10 6 

0.301 

AWCYL-4-1 

[±45 / 0 4 / +45]^ 

7.994 

0.0768 

0.00480 

19.259 10 6 

1.64 ■ 10 6 

0.906 -10 6 

0.299 

AWCYL-5-1 

[+45 / 90 4 / 445] 5 

8.007 

0.0795 

0.00496 

18.615 -10 6 

1.64 ■ 10 6 

0.875 10 6 

0.300 

AWCYL-1 1-1 

[±45/0/90] 2j 

8.0 

0.08 

0.005 

18.5- 10 6 

1.64 • 10 6 

0.87 10 6 

0.3 

AWCYL-92-01 

[±45 / 0 2 ] s 

8.0 

0.04 

0.005 

18.5- 10 6 

1.64 10 6 

0.87 - 10 6 

0.3 

AWCYL-92-02 

[±45/90 2 ] 5 

8.0 

0.04 

0.005 

18.5 10 6 

1.64- 10 6 

0.87 - 10 6 

0.3 

AWCYL-92-03 

[±45/0/90] 5 

8.0 

0.04 

0.005 

18.5 10 6 

1.64 • 10 6 

0.87 - 10 6 

0.3 


For all shells: L = 14.0 in 
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Table 3. Theoretical and experimental buckling loads of axially compressed composite shells. 9 



N c/ [Ib/in] 

K 

y 

Pa 

yP c [ib] 

P a P c [lb] 

P e *p [lb] 

AWCYL-1-1 

-2238.3 

0.3660 

0.471 

0.68 

-19393 

-27999 

-30164 

AWCYL-2-1 

-8910.1 

0.2225 

0.580 

0.68 

-57875 

-67853 

-73975 

AWCYL-3-1 

-8878.7 

0.3980 

0.579 

0.68 

-102905 

-120856 

-147759 

AWCYL-4-1 

-8599.1 

0.3186 

0.575 

0.68 

-79123 

-93572 

-125542 

AWCYL-5-1 

-8892.8 

0.3397 

0.580 

0.68 

-88154 

-103353 

-91662* 

AWCYL-1 1-1 

-8957.4 

0.3983 

0.581 

0.68 

-104182 

-121934 

-152100 

AWCYL-92-1 

-2239.3 

0.2659 

0.471 

0.68 

-14096 

-20351 

-27788 

AWCYL-92-2 

-2239.3 

0.2919 

0.471 

0.68 

-15476 

-22343 

-31924 

AWCYL-92-3 

-2360.4 

0.3561 

0.471 

0.68 

-19901 

-28731 

-34180 


• premature material failure 



Fig. 1 Measured initial shape of a [±45 / 0 / 90] 2s graphite-epoxy shell [9] 
(Radius: 203.2 mm, length: 355.6 mm, wall-thickness: 2.032 mm). 
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refiobMity function R(rho) 


Reliability function for group of 9 NASA composite shells 
Wbor/h 1 fCOs(2*‘pI*x/L) *4- X2*»In (pi^x/L) * cos ( S^y/R) 

FOSM Mont* Carlo 5000 



Fig. 2 Reliability of the group of nine composite shells [9] calculated 
using Koiter’s modified 2-mode imperfection model. 
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